A note about the ground state of the hydrogen molecule 

A.V. TurbineiQ and N.L. GuevarsS 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 

Abstract 

A trial function is presented for the H2 molecule which provides the most accurate (the lowest) 
Bohr-Oppenheimer ground state energy among few-parametric trial functions (with < 14 parame- 
ters). It includes the electronic correlation term in the form ~ exp (7ri2) where 7 is a variational 
parameter. 
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I. INTRODUCTION 



Hydrogen molecule H2 is among the most important chemical objects which appear in Na- 
;ure. Since early days of quantum mechanics after pioneering paper by James and Coolidge 
l| many studies of H2 were carried out (see [2] and references therein). The paper pj] con- 
tained a clear indication that the interelectron correlation must be included explicitly. In 
general, the success of calculations and, in particular, a rate of convergence of a method used 
depends very much on a form of the correlation factor (for a review, see |6i]. Section 

2.2). In particular, recently, it was drawn a conclusion based on an analysis of many atomic 
and molecular systems that the best form of correlation factor is exp(7ri2) comparing to 
linear or the Gaussian in ri2 factors • No clear reason was given so far why it is so. A goal 
of this note is to present a simple, compact, easy-to-handle trial function which leads to the 
most accurate (the lowest) Bohr-Oppenheimer ground state energy among few-parametric 
trial functions (< 14 parameters). The variational energy is calculated numerically using 
a specially designed computer code for multidimensional numerical integration with high 
accuracy. The trial function contains interelectron correlation in the form exp(7ri2). It is 
worth mentioning that long time ago this dependence on ri2 appeared in the variational 



metic fie 
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and, recently, of other 



ll| . A hint why namely this 



trial functions in studies of the H2 molecule in a mag 
two-electron molecular systems in a magnetic field 
ri2-dependence leads to the fast convergent results will be given. 

The Hamiltonian which describes the hydrogen molecule under the assumption that the 
protons are infinitely massive (the Born-Oppenheimer approximation of zero order) can be 
written as follows 

^ =Ep£' - E — + - + I' (1) 

tr rr. ^^''^ ^ 

K=A,B 

where pe = —iVe is the 3-vector of the momentum of the £th electron, the index k runs 
over protons A and B, « is the distance between ith electron and k proton, ri2 is the 
interelectron distance, R is the interproton distance. It is the established fact that the 
ground state of the H2 molecule is the spin-singlet state, symmetric under permutations 
of electron positions as well as proton positions. 



2 



II. VARIATIONAL METHOD 



The variational procedure is used as a method to explore the problem. The reci pe of 
choosing the trial function is based on a physical relevance arguments (see, e.g. 12|). 
In practice, use of such trial functions implies the convergence of a special form of the 
perturbation theory where the variational energy is the sum of the first two terms. Let us 



remind the essentials of this perturbation theory (for details, see 12|). Let us assume that 
our original Hamiltonian has a form Ti. = — A + V. As a first step we choose a trial function 
^(triai) is normalized to one. Then we find a potential for which our trial function 

^(triai) exact eigenfunction Vtriai = ^^^^0' with the energy Etriai = 0. In a pure formal 

way we can construct a Hamiltonian Htriai = —A + Vtriai such that Ti.triai'^^^"''^'' = 0. It can 
be easily shown that the variational energy 

is nothing but the first two terms in the perturbation theory where the unperturbed problem 
is given by Titriai and the perturbation is the deviation of the original potential V from the 
trial potential Vtriau namely, Vperturbation = V — Vtriai- Eventually, we arrive at the formula 

Euar Etriai ~l~ EiiVpg^rturbation) j (2) 

here Ei(yperturbation) =< '^^^'''''^^ \Vperturbation\i^^^''"^^^ > IS the first energy correction in the 
perturbation theory, where unperturbed potential is Vtriai- It is worth noting that if the 
trial function is the Hartree-Fock function the resulting perturbation theory is nothing but 
the Moeller-Plesset perturbation theory (see, e.g. 13 1, Section 15.18) [2^. 

One of the criteria of convergence of the perturbation theory in Vperturbation = V — Vtriai 
is a requirement that the ratio | ^p^rturbauon | gj-^Q^j^j j^q^ grow when r tends to infinity in any 
direction. If this ratio is bounded by a constant it should be less than one. In fact, it is 
a condition that the perturbation potential is subordinate with respect to the unperturbed 
potential. A value of this constant controls the rate of convergence - a smaller value of this 



constant leads to faster convergence 15|. Hence, the above condition gives a importance to 
the large-range behavior of the trial functions. In the physics language the above requirement 
means that the phenomenon of the Dyson's instability should not occur (for a discussion 



see 



12 1) 21 1 . Among three factors which are mentioned in literature (see jfi]): the linear 
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in ri2, exponential exp(7ri2) and expf— arjg)) the only factor exp(7ri2) fulfills the above 
requirement. It was demonstrated in |3| that a superposition of the Coulomb functions with 
exponentially correlated function exp(7ri2) (see below, eq.(4)) leads to faster convergence 
than others. Perhaps, it is worth mentioning that for the case of Gaussian factor the above- 
defined constant is equal to one exactly. In concrete, by following the above procedure and 
a requirement of the convergence of the perturbation theory we choose the trial function for 
the ground state in a form 

^itr^al) ^ ^^^^ ^ ^^^^ ^ ^^^^ (3) 

where 

-02 = (1 + P-^2)e~"'^^'^^"^~^^^^'^~"^^^^^~^^'^"^'^^'^'^^^'^ (5) 

— (•]_ _|_ p_^^''jQ-a7iriA+riB)-as(r2A+r2B)+-y3ri2 ^ ^g-J 

The P12 is the operator which interchanges electrons (1^2) and Pab is the operator 
which interchanges the two nuclei A ^ B. It is easy to check that the functions (5)-(6) are 
symmetric with respect to the interchange A ^ B. The variational parameters consist of 
non- linear parameters ai_8, 71-3 which characterize (anti) screening of the Coulomb charges 
and linear parameters Ai^^. If the internuclear distance R is taken into account the trial 
function ([3]) depends on 14 parameters 



22 |. It is worth mentioning that ([5]) is a degeneration 
of (jlj) when ai = 04, 0^2 = ctz and ([6]) is another degeneration of (jl]) when ai = 02, "3 = 04. 
In a certain way, the function ([5]) mimics the interaction of two hydrogen atoms H + H, 
while the function ([6]) mimics the interaction H2 + e. Eventually, the function (jlj) can be 
treated as a non-linear interpolation between ([5]) and ([6]). Those functions look analogous 
to the Hund-Mulliken, Heitler-London and Guillemin-Zener functions, respectively. 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Multidimensional integration was carried out numerically using a " state-of-the-art" dynam- 
ical partitioning procedure: a domain of integration was divided into subdomains following 
an integrand profile, in particular, the domains with sharp changes of the integrand were 
separated out. Then each sub domain was integrated separately with controlled accuracy 
(for details, see e.g. [3]). A realization of the routine requires a lot of attention and care. 
During minimization process a partitioning was permanently controlled and adjusted. Nu- 
merical integration was done with a relative accuracy of ~ 10^^ — 10^ ' by use of the adaptive 



DOIFCF routine from NAG-LIB. Computations were performed on a dual DELL PC with 
two Xeon processors of 2.8 GHz each. 



III. RESULTS 

Present results for the ground state of the H2 molecule and their comparison with results 
of previous studies are displayed in Table [H The Bohr-Oppenheimer ground state energy 
obtained using the function (3)-(6) is the most accurate (the lowest) among those obtained 
with other trial functions with < 14 parameters. A reasonable agreement for expectation 
values is also observed, except for < 3zf — rf > related to the quadrupole moment. To 
present authors it seems evident that this expectation value should be studied separately 



(see also 17|). It is not surprising that the obtained value of the cusp parameter in ri2 is 
equal to 0.4 unlike to the exact value 0.5 (see footnote [19]). Variational parameters of the 
trial function ([3]) are shown in Table HTl It is worth emphasizing that a numerical calculations 
are very difficult and can easily lead to a loss of accuracy. In jsl a similar function (3) but 
with all three components of the form (4) containing 18 variational parameters was studied 
using variational Monte-Carlo technique. A comparison of our results with less parameters 
with 3|] (see Table I) indicates that we obtain a lower total energy of the order 5 x lO^'* Ry. 

IV. CONCLUSION 

We presented a simple and compact few-parametric trial function which provides the most 
accurate Bohr-Oppenheimer energy for H2 molecule among those based on few-parametric 
(< 14) trial functions. Emerging five- dimensional integrals were effectively calculated using 
fast state-of-the-art integration routine which admits parallelization. The trial function 
(j3]) can be easily generalized by adding other physically-natural degenerations of (j4]) than 
([5]), ([6]). One of them appears when all a-parameters in (jl]) are equal. It should be dominant 
in a domain of small interproton distances. It seems natural to assume that taking linear 
superpositions of the functions (jlj) we end up with fast convergent procedure (see [3|]). The 
function ([3D can be easily as modified for a study of spin-triplet states and as well as the 
states of the lowest energy with non-vanishing magnetic quantum numbers. A generalization 
to more-than-two electron systems is straightforward. 
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TABLE I: Total energy Ex in Ry and expectation values in a.u. of the hydrogen molecule H2 for 
the ground state. ri,zi are distances from 1st electron to the mid-point between protons. Some 
data are rounded. 



Et (Ry) 


< ^12 > 


< rf > 


R 


< 3zj - rj 


> 


Refs. 


-2.34697 












[1] 


-2.34778 ^ 












[18] 


-2.34787 " 












[3] 


-2.348382 




2.5347 




0.5227 




[17] 


-2.348393 ^ 


0.5874 


2.5487 


2.2133 


0.4847 




present 


-2.34872 




2.5426 




0.5142 




m 


-2.34888 3 


0.5874 




2.2127 






m 


-2.34895 ^ 












m 



(the BO energy with 14 variational parameters) 



"from t(" 

^from T|8| (based on use of > 200 non-spherical Gaussian orbitals) 
"Exponentially correlated functions) 



'^from [3| (A^ = S 
''from Table III 
•^Present calcula' 



17| (the BO energy with 14 variational parameters) 



^from 
''from 



TUti (14 variational parameters) 



■''from Table lIlrTF| (the BO energy with 28 variational parameters ) 



17| (the BO ground state energy with 40 variational parameters) 



ble II 

(7034 James-Coolidge type terms, the record calculation at present, the number in Table is 



rounded)) 
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TABLE 11: Parameters of the trial function ([3]). 



R 


1 40053 




1. 




720fi74Q8fi 




521577488 


Ol'X 


0.130799743 




1.30816746 


n 


0655006315 




-1 15105579 




0.604583808 




658402827 


72 


-0.349101361 


^3 


0.256342676 


a^ 


0.968330781 


as 


0.229153253 


73 


-0.354509413 
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